ESC  1  l<-2(M)7-()(.7 


Technical  Report 
1135 


I 


I 


r 


I 


‘‘Phase-Enhanced”  3D  Snapshot  ISAR 
Imaging  and  Interferometric  SAR 


J/r.  Mayhan 


28  Doccmhcr  2009 


Lincoln  Laboratory 

MASSACin  SETTS  LNSn  iT  TK  OF  TFCllNOFOGY 

Lh.xi \(;r()\,  M \ ns u  / / 1  shiis 


Prepared  for  tlie  Missile  Defense  Agency  iiinler  Air  Force  (^onlrac!  FAB72l-()!>-C-(HM)2. 


Approved  for  pnl)lic  l  elease:  distri!)iition  is  iinliniited. 


20100105288 


V 


This  report  is  basod  on  studies  performed  at  Lincoln  Laboratory,  a  center 
for  research  operated  by  Massachusetts  Institute  of  Technolofy.  This  work 
was  sponsored  by  the  Missile  Defense  Agency,  DV,  under  Air  Force  G>n tract 
FA8721-05-C-0002. 

This  report  may  be  reproduced  to  satisfy  needs  of  U.S.  Government 
agencies. 


The  MDA  Public  Affairs  Office  has  reviewed  this  report, 
and  it  is  releasable  to  the  National  Technical  Information 
Service,  where  it  will  be  available  to  the  general  public, 
including  foreign  nationals. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


FOR  THE  COMMANDER 


Gary ' 

AdminiVfrative  Contracting  Officer 
Plans  and  Programs  Directorate 
Contracted  Support  Management 


Non-UrKOln  Recipients 

PLEASE  DO  NOT  RETURN 

Permission  has  been  given  to  destroy  this 
document  when  it  is  no  longer  needed. 


i\1assa<‘hiis<‘lls  Insliliil<‘  of  "Pechnolojiv 
Lin<‘olti  Lahoralory 


“■Phas<‘-f]nhan<*<*(l’' 31)  Snapshot  ISAP  linaj^inj^  and  lnl<‘rforoni<‘lrir  SAI? 


J/r.  Mdvlidn 
(}roiif)  32 


Trcliiiical  Hcporl  1  133 


2«  l)<H  <‘ml)er  2009 


Approvt'cl  for  piiblir  rcirase;  <listril>iitioii  is 


l^<*xin'^lon 


Massarhiisolls 


This  page  intentionally  left  blank. 


ABSTRACT 


In  [1],  we  presented  a  novel  fomiiilalion  for  the  generation  of  three-dimensional 
(3D)  inverse  synthetie  aperture  radar  (ISAR)  images  based  on  reeent  developments  in  high 
resolution  speetral  estimation  theory'.  Beeause  this  teehnique  requires  only  '“’snapshots’'  of 
data,  where  a  '"'snapshot’'  is  defined  as  a  bloek  of  Nyquist  sampled  freqiieney-time  pulses 
of  data,  we  refer  to  it  as  '‘3D  snapshot  imaging.'’  Coneomitant  with  these  results,  reeent 
advanees  in  interferometrie  SAR  imaging  have  demonstrated  the  use  of  two-dimensional 
(2D)  range-Doppler  image  phase  information  to  extraet  '‘out-of-plane"  height  information 
to  obtain  3D  images  of  ground  seenes.  In  this  ease,  a  unique  sampling  grid  is  generated 
that  allows  ‘'overlaying’'  of  nearly  identieal  2D  range-Doppler  images  and  uses  phase 
diff'erenees  between  these  images  to  estimate  the  “out-of-plane"  height  infonnation,  from 
whieh  a  3D  image  is  developed. 

In  this  report,  we  develop  a  framework  eonneeting  these  two  teehniques, 
partieularly  applieable  to  forming  3D  images  of  target  types  typieally  dominated  by 
smaller  numbers  (<20)  of  seattering  eenters,  and  eharaeteri/ed  by  detenninistic 
exoatmospherie  motion  having  torque-free  Euler  dynamic  spin  and  preeession.  For  the 
application  considered  here,  the  3D  image  is  formed  using  the  local i/ed  mierodynamic 
changes  in  look  angles  to  the  target  eharaeteri/ed  by  Euler's  6  Degrees-of-Freedom 
(6DOF)  dynamical  equations  of  motion.  Fhis  type  of  motion  allows  the  formation  of 
short-time-interval  range-Doppler  images  at  nearly  constant  look  angle  Kappa  (K)  to  the 
time  invariant  angular  momentum  vector  J.  This  short-time  processing  is  contrasted  with 
longer-time  ISAR  imaging  of  satellites,  where  changes  in  look  angle  to  the  target  are 
obtained  from  changes  in  the  long-term  trajectory'  of  the  satellite. 

Applications  of  the  development  are  illustrated  using  simulation  data;  speeifieally, 
the  w  ideband,  monostatie  field  is  simulated  using  a  simple  point  scatter  model  of  a  generic 
eone-like  target,  and  the  phase-enhanced  3D  image  is  generated  for  differing  cases. 
Extensions  of  the  teehnique  to  bistaties  are  also  discussed. 
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1.  INTRODUCTION 


Three-dimensional  (3D)  radar  imaging  of  inverse  synthetie  aperture  radar  (ISAR)  data  presents  a 
ehallenge  of  eonsiderable  importance  to  numerous  application  areas.  The  development  and  study  of  3D 
ISAR  imaging  and  associated  etTects  has  been  carried  out  by  various  researchers  [2-4].  A  general 
formulation  of  this  3D  imaging  problem  is  presented  in  [2|.  The  essence  of  that  fonniilation  is  to 
eharaeteri/e  the  target  by  setting  up  a  3D  grid  registered  in  a  coordinate  system  centered  within  the  target, 
and  assoeiating  two-dimensional  (2D)  range-DoppIcr  image  planar  cuts  through  the  target  with  points  on 
this  grid  using  the  (assumed)  known  aspect  angle  to  the  target.  When  the  radar  samples  a  planar  cut 
through  the  target,  the  resulting  image  will  necessarily  be  2D.  When  the  radar  densely  samples  a  3D  solid 
angle  of  target  aspects,  a  3D  image  of  these  seatterers  can  be  developed. 

In  a  previous  paper  [I  ],  we  presented  a  new  formulation  for  the  generation  of  3D  ISAR  images  based 
on  recent  developments  in  high  resolution  spectral  estimation  theory.  Because  this  technique  requires  only 
high  resolution  ‘'snapshots’'  of  data,  we  refer  to  it  as  ‘'3D  snapshot  imaging."  Because  of  the  snapshot 
nature  of  the  technique,  it  is  particularly  applicable  to  3D  imaging  of  sectors  of  sparse  angle  data. 
Concomitant  with  these  results,  3D  SAR  imaging  of  ground  refleeting  surfaces  has  shown  much  promise 
for  synthetic  aperture  radar  processing.  A  key  factor  in  achieving  such  3D  images  is  precise  control  over 
the  grid  (motion)  used  by  the  radar  sensor.  A  particular  advaneement  in  this  research  area  is  the  technique 
referred  to  as  interferometric  SAR  (IF-SAR).  For  this  ease,  a  unique  sampling  grid  is  generated  that  allows 
“overlaying"  of  nearly  identical  2D  range-DoppIcr  images  and  uses  phase  difTerenccs  between  these 
images  to  estimate  the  “out-of-plane"  height  information,  from  which  a  3D  image  is  developed. 

In  this  report,  we  develop  a  framework  eonneeting  these  two  techniques,  particularly  applicable  to 
forming  3D  images  of  target  types  typically  characterized  by  smaller  numbers  (5  20)  of  dominant 
scattering  centers,  and  characterized  by  deterministic  exoatmospherie  motion  having  Buler  dynamic  spin 
and  precession.  For  the  application  considered  here,  the  3D  image  is  fonned  using  the  localized 
microdynamie  changes  in  look  angles  to  the  target  characterized  by  Euler’s  6  Degrees-of- Freedom  (6DOF) 
dynamical  equations  of  motion.  This  type  of  motion  allows  the  formation  of  short-time-interval  range- 
Doppler  images  at  nearly  constant  look  angle  Kappa  (K)  to  the  time  invariant  angular  momentum  vector  J. 
This  short-time  processing  is  contrasted  with  longer-time  ISAR  imaging  of  satellites,  where  changes  in 
look  angle  to  the  target  are  obtained  from  changes  in  the  long-term  trajectory  of  the  satellite. 

In  essence  a  SAR  interferometric  radar  measures  height  by  comparing  the  phase  of  corresponding 
range- Doppler  image  planes,  calibrated  for  sensor  motion,  to  estimate  the  ‘'projection  height"  or  3D 
information  contained  in  the  two  nearly  identical  in  amplitude  data  sets.  References  [5-7]  provide  a  good 
discussion  of  the  basic  principles  and  associated  radar  calibration  requirements.  Various  schemes  for 
eolleeting  the  SAR  data  arc  discussed,  each  focused  on  the  peculiar  hardware  calibration  requirements 
inherent  in  applying  these  techniques.  For  the  ISAR  sensor  development  presented  here,  two  cases  are 
considered:  monostatie  (single  sensor)  and  bistatic  (two  sensors).  When  applied  to  a  single  sensor,  these 
calibration  problems  arc  essentially  eliminated  and  replaced  by  the  need  for  an  accurate  motion  solution 
and  apply  to  a  specific  type  of  target  motion.  For  the  bistatie  case,  a  different  set  of  calibration 
requirements  apply,  as  discussed  in  subsequent  text. 


For  a  single  monostalic  sensor,  the  essence  of  the  applicability  of  these  IF-ISAR  techniques  to  Euler- 
type  motion  is  the  observation  that,  for  the  case  when  the  spin  frequency,  f^,  is  much  greater  than  the 
precession  frequency,  fp,  and  if  the  precession  cone  angle  is  not  too  large,  range-Doppler  images  at  roll 
intervals  (j)  and  0  +  271  look  nearly  identical,  making  the  range-Doppler  scatterer  correlation  problem 
relatively  simple.  Once  individual  scattering  centers  are  correlated  between  images,  the  phase  difference 
between  images  can  readily  be  determined  and  used  to  augment  the  basic  snapshot  imaging  equations 
developed  in  [1].  (Note:  In  theory,  the  phase  difference  between  scattering  locations  can  be  used  at  any 
point  along  the  dynamical  motion.  However,  the  practical  problem  of  tracking  and  correlating  specific 
scattcrers  over  extended  time  intervals  generally  prohibits  this.)  We  show  that  the  utility  of  the  technique 
degrades  gracefully,  i.c.,  when  the  precession  cone  angle  becomes  large,  and  the  increment  in  aspect  A6 
occurring  at  the  roll  angles  0  and  0  +  271  is  greater  than  allowed  by  Nyquist  sampling  (for  which  the 
extracted  phase  difference  between  images  for  some  scatterers  exceeds  271,  and  therefore  becomes 
ambiguous),  the  standard  snapshot  equations,  without  phase,  remain  valid;  yet  one  still  has  full  advantage 
of  the  simpler  track  correlation  problem. 
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2.  EULER  MOTION  AND  SAMPLING  CONSIDERATIONS 


We  assume  a  target  eentered  eoordinate  system  as  illustrated  in  Figure  1 . 


Target 


To  radar 

Radar  line  of  sight 
changes  with  time 


y 


Figure  I.  Motiou  and  target  scattering  reference  frames. 


Standard  spherical  coordinates  6,  ({)  eharaeteri/e  the  look  angle  to  the  radac  which  is  determined  by 
the  unit  vector  k,  where  the  notation  denotes  a  unit  vector.  In  general,  6(t)  and  (t)(t)  are  complex 
functions  characterized  by  Euler's  dy  namical  equations  of  motion.  The  target  coordinates  (and,  hence,  the 
image  reference  frame)  are  defined  relative  to  the  origin  of  this  system,  and  an  arbitrary'  scattering  point  on 
the  target  is  denoted  as  (Xp,  y^,  z^).  For  a  wideband  radar  operating  in  the  higher  frequency  bands,  the  net 
scattering  from  the  target  can  be  decomposed  into  scattering  from  a  number  of  discrete  scattering  centers, 
denoted  by  (Xp,  y^,  Zp),  n=l,...,  N.  Assume  we  can  isolate  a  specific  seatterer  rp  =  (Xp,  yp,  Zp)  on  the  target. 
In  a  target-centered  coordinate  system,  this  position  is  fixed,  and  the  look  angle  to  the  radar,  and  its  time 
derivative,  are  characterized  by 


k  •  r 


n 


=  R.. 


(1) 


k-r_ 


(2) 


where  we  define  (R„,  R„)  to  be  the  range,  range  rate  radar  observables  eorresponding  to  (Xp,  z„).  and  k 
and  k  are  funetionsof  0,d.(j),(j)  given  by 
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k  =  sin  6  cos  (j)  x  +  sin  6  sin  (f)  y  +  cos  6  2 


(3) 


k  =  66  +  ^  sin  6  (p 


(4) 


where 


6  =  cos  6  cos  (p  X  +  cos  6  sin  (p  y  -  sin  6  2 
(p  =  -sin(px  +  cos  (p  y 


Eqs.  (I)  and  (2)  form  a  nonlinear  set  of  equations  eharaeterizing  the  behavior  of  0(t)  and  (l)(t)  and  eaeh 
seattering  eenter  (x^,  y^,  7,„).  We  impose  physieal  constraints  on  the  motion  0(t)  and  (l)(t)  by  constraining 
0(t)  and  <|)(t)  to  that  body  motion  following  from  Euler’s  dynamical  equations  in  a  torque-free  environment. 
To  this  end,  we  require  a  closed  form  transformation,  E,  of  the  type 


F[K,6p,(Pp,y/,f,,fp,a,,ap,t]^6{t)and  (p{t) 


(5) 


where 


0p 

¥ 

fs 

fn 


localized  aspect  angle  from  the  radar  to  the  target  angular  momentum  vector 

precession  cone  angle  relative  to  target  angular  momentum  veetor 

precession  rotation  angle 

spin  angle 

spin  frequency 

precession  frequency 

spin  phase 

precession  phase 


Bq.  (5)  imposes  physieal  constraints  on  the  motion  valid  for  torque-free  exoatmospherie-type  targets. 
The  elosed-form  relationship  defined  by  (5)  is  readily  determined  using  the  coordinate  transformations 
developed  in  [8]  and  is  summarized  below. 


The  aspect  angle  0(t)  is  dependent  on  the  precession  parameters  according  to 


0(t)  ^cos~^  ^cosKcosO p  ^siriKsinO pSinipp  ^ 


(6) 
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The  roll  angle  ^(t )  is  dependent  on  spin  and  precession  variables  according  to 


0{t)  =  -\-y/  -tan  ^ 


cos  Ksin  6^^  -  sin  Kcos  0^^  sin0^, 
sin  K’cos 


(7) 


The  elosed’form  analytic  expressions  defined  by  (6)  and  (7)  arc  cxtremcl>  useful  in  understanding 
the  dependence  of  (0,  d))  on  any  specific  type  of  Euler  motion. 

A  special  case  of  interest  occurs  for  0p  small.  In  this  case,  the  expressions  0(t),  (t)(t)  reduced  to  the 
simpler  forms  [carrying  out  the  small  angle  expansions  to  order  (0p)] 


Oft  )  =  K-0pSin0^  ' 


(8) 


0(t  )  =  271  f^t  -  dp  ctnKcos  0p  ^ 


(9) 


where 


(Pp=27ifpt  +  ap. 


(10) 


and  we  define  the  ‘"first  order"  roll  frequency  f^^  =  f^  +  fp  =  where  is  the  first  order  roll  period, 

fhe  roll  period,  f^ ,  is  in  fact  most  directly  inferred  from  radar  observables.  Note  that  as  defined  above 
is  a  first-order  approximation  to  the  instantaneous  roll  frequency  T(t)  =  2k/^  . 

Figure  2  illustrates  the  temporal  variation  of  0(t),  (t)(t)  and  the  corresponding  vector  k(t)  for  the  ease 
where  0p  =  2.7^  f^  ~7  *  fp,  and  Tp  =  4  seconds.  Note  the  helical  nature  of  k(t),  plotted  over  the  time 
interval  Tp/2,  characteristic  of  spinning,  prccessing  motion.  We  can  use  the  variation  of  k(t)  in  Figure  2  to 
contrast  the  sampling  grid  that  occurs  using  traditional  IF-SAR  techniques.  This  contrast  is  illustrated  in 
I'igure  3. 
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Figure  2.  Euler  motion  (0,  (pcmd  k  time)  for  small  precession  and  f  » fp. 


Time  Variation  Over  Interval  Tp/2 


Sampling  Grid  for 
Conventional  IF  -  SAR  Imaging 


0) 


Sampling  Grid  for 
IF  -  ISAR  Euler  Motion 
fs  >  >  fp 
0  small 

p 


0  In  Phi 


Figure  3.  IF- ISAR  sampling  grids. 

Note  that  for  conventional  IF-SAR,  the  aspect  0  is  constant  over  the  principle  (0,  2k)  0  intervals 
shown,  whereas  for  spinning,  prccessing  motion,  0(t)  is  constantly  changing  over  the  sampling  interval.^ 
For  either  case,  extracting  an  unambiguous  phase  difference  between  regions  ([)  and  0  +  271  requires  the 
sampling  interval  A0|vf  be  less  than  Nyquist,  i.e.. 


*  The  sampling  grid  illustrated  for  IF-ISAR  for  spin  and  precession  motion  is  somewhat  simplistic,  for  illustration 
purposes.  See  Figure  10  for  a  more  typical  ease. 
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<Xq/  (2*  Projeeled  largel  Length), 


(M) 


where  denotes  the  wavelength  at  hard  eentcr  1'=  Tq. 

As  mentioned  earlier,  the  utility  of  using  3D  snapshot  imaging  is  the  ability  to  generate  3D  target 
images  over  only  a  sparsc-angle  observation  spaee.  In  Figure  4,  we  eompare  several  eases  illustrating  the 
physieal  oeeurrenee  of  sparse-angle  sampling  eomparcd  to  the  full  4*7r  steradian  sampling  spaee.  The 
sparseness  of  spin,  precession  sampling  is  clearly  evident.  For  comparison,  RCS  measurement  data 
collected  on  a  typical  radar  cross  section  (RCS)  range  most  clearly  correlates  to  the  special  cases  of  "pure 
tumble''  ((t)  =  constant)  and  ‘'pure  spin"  (0  =  constant).  The  sampling  grid  associated  with  phase-enhanced 
IF-ISAR  for  the  special  case  fs»fp  considered  in  this  report  is  also  shown  for  comparison. 
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figure  4.  Types  of  sparse-angle  sampling. 


Figure  5  contrasts  the  sampling  grid  behavior  using  phase-enhanced  snapshot  processing  vs. 
conventional  snapshot  processing.  Using  conventional  snapshot  processing,  vve  sec  that,  for  f^  »  fp,  the 
spectral  sequence  variable  R  is  dominated  by  (j)  information  over  a  single  roll  segment,  whereas  phase- 
enhanced  snapshot  processing  provides  both  0  and  0  Doppler  information  by  including  the  0  samples 
between  fixed  roll  cuts.  We  show  later  that  this  results  in  considerably  greater  linear  independence  of  the 
phase-enhanced  snapshot  equations  vs.  conventional  snapshot  imaging. 
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Conventional  Phased  Enhanced 

Snapshot  Imaging  Grid  Snapshot  Imaging  Grid 


Figure  5.  Snapshot  sampling:  conventional  V5.  phase-enhanced. 
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3.  DEVELOPMENT  (SINGLE  MONOSTATIC  SENSOR) 


Consider  the  geometrv  illustrated  in  I  igure  I.  Using  the  Cieomelrieal  Theorv  of  Diffraction  ((i TD) 
[9],  applieable  to  the  high  freqiieney  region  under  eonsideration,  monoslalie  baekseatter  from  the  target 
ean  be  written  in  the  form 


k,  .  r„ 


(12) 


where  the  set  {r^},  N  denotes  the  N  dominant  scattering  centers  on  the  target,  each  located  at 

(Xp,yp,/p),  having  complex  amplitude  Dn,  and  denotes  the  radar  line-of-sight  (LOS)  unit  vector. 
Typically,  Dn  is  a  relatively  slow  var>'ing  function  of  (f,  0,  (}))  and  eharaeteri/.es  the  GTD  difTraelion 
coefficient  of  the  nth  scattering  center.  It  is  essentially  constant  over  the  dala-bloek  under  consideration 
when  using  sequential  block  processing.  Denote  k  =  k ,  at  lime  I  =  t|,  and  k  =  k  ^  at  t  -  t2  where  I2  >  t|. 
Now  assume  that  at  the  later  lime,  t  =  I2,  the  scattered  field  is  sampled  in  a  direction  k  .  .  Furthermore,  if 
we  assume  k  ^  is  "close"'  to  k ,  ,  then  we  ean  express 


k  ,  =  k  , + Ak  • 


(13) 


Using  (13)  in  the  general  expression  for  the  seallered  field,  the  scattered  field  1:  ^  at  t  =  t2  can  be 
expressed  in  the  form 


I'-2=S|D1  c 


k  .  r 
J  -i  ‘11 


(14) 


Now  assume  that,  since  k .  is  assumed  close  to  k,  .  the  2D  range-Doppler  images  extracted  from 
each  2D  data  block  centered  at  times  t  +  t)  and  t  =  I2,  respectively,  are  nearly  identical,  so  that  the 
individual  scattering  centers  in  each  image  ean  readily  be  identified  and  correlated.  By  extracting  the  phase 
of  each  scattering  center  from  the  range-Doppler  image,  Eqs.  (12)  and  (14)  ean  be  compared  to  obtain 


Phaseilj^) 


/V7a.se(/|„)  =  Phuse{D\)-Phase{D^^)  +  Y(^2-k\)*r„ 


(15) 
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where  Ij^  denotes  the  complex  amplitude  of  the  image  pixel  of  the  nth  scattering  center.  If  k2-k\  is  close 
to  Nyquist  sampling,  then  D  D  ^  so  that  ( 1 5)  can  be  written  in  the  form 


4tc  ^  ^ 

Kn  =  ^„{Phase), 


(16) 


where  we  define: 


/S.^{Phase)  =  Phase{l2n)-  Phase{I^^) 


(17) 


where  Ijp  and  l2n  represent  the  complex  amplitude  defined  over  each  respective  range-Doppler  image 
plane. 


Eq.  (16)  will  be  used  in  the  next  section  to  augment  the  standard  set  of  3D  snapshot  imaging 
equations. 

It  should  be  noted  that  Eqs.  (15-17)  as  derived  are  equivalent  to  the  standard  characterization  of  IF- 
SAR  processing.  They  are  repeated  in  this  manuscript  only  to  emphasize  their  coupling  to  the  3D  snapshot 
imaging  equations. 
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4.  THE  GENERALIZED  ^THASE-ENHANCED’’  SNAPSHOT 
IMAGING  FRAMEWORK 


In  this  section,  we  first  review  the  methodology  of  generating  a  3D  snapshot  image  described  in  [  1 1 
and  then  augment  this  basic  equation  set  with  the  phase-enhanced  equations  developed  in  Section  3. 

Consider  First  the  process  of  sequentially  extracting  the  range,  range-rate  information  from  a  2D 
radar  data  block.  Typically,  a  wideband  radar  transmits  a  sequence  of  pulses,  each  of  which  interrogates  the 
target  over  a  given  frequency  band  fQ  -  BW/2  to  f()  +  BW/2,  where  f()  is  the  carrier  frequency  and  BW 
the  bandw  idth  of  the  radar.  Haeh  pulse  in  the  time  sequence  is  spaced  at  time  interv  al  At,  where  the  PRF  of 
the  radar  is  given  by  PRF  =  1/At.  Denote  the  signal  received  (after  matched  filtering)  by  the  radar  on  each 
pulse  as  Cq  (t),  where  the  index  "‘q'"  denotes  the  pulse  label.  Denote  the  unweighted  Fourier  transform  of 
Cq  (t)  as  F^q  (CO),  where  O)  =  271  f  is  a  eontiniious  frequency  variable  and  f  ranges  over  the  bandwidth  (-F3\\7 
2  to  BW/2).  Assume  f  is  sampled  at  frequency  intervals  f  “  k  Af  Then  the  data  matrix  F]q  (k)  “  F  (k,  q) 
represents  a  2D  (frequency,  pulse)  representation  of  the  field  scattered  by  the  target.  A  conventional  linear 
image  could  then  be  obtained  using  a  2D  Fourier  transform  on  E  (k,  q).  Vor  our  purposes,  it  is  more  useful 
to  consider  the  raw  data  block  F  (k,  q)  and  optimally  extract  target  information  from  this  block  using 
recently  developed  2D  spectral  estimation  techniques  [10,1 1,12].  As  time  "'marches  on,"'  sequences  of  data 
blocks  each  incremented  by  At  can  be  processed.  We  refer  to  this  as  sequential  block  processing.  Fhe 
overall  process  is  illustrated  in  Figure  6  using  simulated  data  from  a  cone- like  target  having  extended  fins 
(see  Figure  8). 


Single  Block  Spectral  Estimate 


1 


N=  1 1  Scatterers 


-1,5 


-3  -2  0  1  2  3 

Range  Rate 


Figure  6.  Compoueut  extraction  using  block  processing  vs.  conventional  Fourier  image. 


Two-dimensional  blocks  of  radar  data  are  processed  sequentially  to  develop  a  range,  range-rate 

spectral  sequence  =  /y . as  illustrated  in  Figure  6.  Typically,  an  R,R  correlator/smoother  is 

implemented  to  track  these  poles  in  R^R  space  by  sequentially  processing  a  block  of  pulses  through  the 
measurement  data  set.  Eqs.  (1,2)  are  then  applied  sequentially  to  a  number  of  snapshots  to  develop  an 
image  of  the  scattering  points  .  7  b  illustrate  this  process  mathematically,  express  the 

transformation  defined  by  Eqs.  (1)  and  (2)  in  matrix  form: 


T 

=^// 


n 

y,, 

A. 

Jn. 

(18) 


where  T_  is  a  2  x  3  transformation  matrix  given  by 


where  k  and  k  are  expressed  as  row  vectors  and  (18)  defines  the  mapping  at  t  tqj.  We  can  accumulate 
these  equations  for  two  snapshots  at  tqj  and  tq2  in  the  form 


—  - 

X 

r 

='// 

T 

yn 

y 

V 
/ 

V 


R  ^ 

R 

R  ^ 
R 


nj 


(19) 


As  additional  snapshots  are  obtained,  these  equations  are  added  to  (19),  which  is  then  solved  in  a 
least  squares  sense  for  the  scattering  coordinates  (x^ » ^  )  • 

The  rank  of  T_  and  hence  the  solutions  for  clearly  depend  on  numerous  variables. 

Some  of  these  arc  the  observed  motion  of  the  target,  the  number  of  snapshots,  the  time  spacing  between 
snapshots,  and  the  signal-to-noise  ratio  present  in  the  data.  A  detailed  analysis  of  these  variables  relative  to 
their  effect  on  the  resultant  image  has  not  been  carried  out.  In  applications  to  actual  field  data  not  presented 
here,  we  demonstrated  good  results  by  implementing  some  smoothing  of  the  range,  range-rate  data  before 
applying  (19). 
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Consider  now  augmenting  the  basic  snapshot  equation  set  with  using  the  phasc-extraeted 
information  defined  b>  Eq.  (16).  Note  that  by  examining  Figure  5,  R  information  is  extracted  from  data 
along  the  direction  of  motion.  For  example,  for  spin-only,  at  constant  aspect  6^^,  0-0 ,  and  using  Fq.  (4) 


<P  sin . 


i.c..  R  is  characterized  solely  by  ‘P  for  this  case.  This  corresponds  to  the  loeali/ed  spin-dominated  motion 
illustrated  in  Figure  5  for  conventional  snapshot  processing  when  f^  »  fp. 

The  sampling  grid  for  phased-cnhanced  snapshot  imaging  illustrated  in  Figure  5  provides  additional 
independent  information  by  augmenting  Eq.  ( 1 8)  with  data  blocks  associated  with  changes  in  0,  i.c.,  6  =  6q 
+  A6  when  0  =  (|)  +  2n.  llq.  (18)  is  then  augmented  by  the  additional  Fq.  (16)  developed  in  Section  3: 


^(k^-ki)»r„=A„  (phase).  (20) 

where  the  subscript  ""2''  is  defined  to  be  at  sample  time  corresponding  to  roll  angle  04-  2/1  and  (phase) 
is  the  phase  difference  for  the  scatter  at  )  at  each  respective  roll  angle.  Clearly  for  A  (phase)  to 

be  nonambiguous,  the  change  A0  in  aspect  cannot  exceed  Nyquist. 

Note,  in  principle,  augmenting  Eq.  (18)  with  (20)  provides  three  independent  equations  to  directly 
solve  for  the  three  unknowns  ).  Fhat  F^q.  (20)  provides  independent  information  can  be  readily 

seen  by  evaluating  -k^)  directly  assuming  0  =  0^+  A0ai  0  =  0  +  2;r  .  It  can  readily  be  shown  that 


471  ''  ''  471 

—  (k,-k,).  I,  =  yAOO 


The  three  resulting  equations  arc  approximately  given  by 


(21) 


sin  6  cos  +  sind  sin0y\^  +ca.v6?z,,  = 


(22) 


0sin00*r^  = 


(23) 
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(24) 


—  A0e*r„=A„  {phase)  ■ 


Eqs.  (22),  (23),  and  (24)  have  an  interesting  graphical  interpretation  which  illustrates  the  basic 
orthogonality  (and  consequent  linear  independence)  of  the  equation  set,  as  illustrated.  The  localized 
orthogonal  system  is  illustrated  in  Figure  7.  Recall  that  k  is  a  unit  vector  in  the  direction  of  the  radar  TOS. 
Since  k  is  perpendicular  to  the  three  unit  vectors  (A:,^,  0)  form  a  local  (i.e.,  a  function  of  6,(l>) 
orthogonal  set  of  unit  vectors,  defining  a  local  3D  axis  set.  Eqs.  (22-24)  then  can  be  written  in  the  form 

k*Ln=Rn 

^•r„=R„/{(psin0]  (25) 

-  A 

9^r_^  -  —  •  A^{phase)/ A0^ 

471 

which  clearly  delineate  the  projections  of  the  scattering  center  location  vector  r^^  onto  the  axes  of  this 
localized  orthogonal  frame. 


Figure  7.  Localized  orthogonal  coordinate  system  for  3D  phase-enhanced  precessing. 
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Although  Eq.  (25)  is  directly  invertible  to  solve  for  )  using  a  single  "'snapshot,"  it  is  still 

useful  to  use  a  cumulativ  e  approach  that  prov  ides  smoothing  in  the  presence  of  noisy  data.  To  this  end  we 
rewrite  Eq.  (25)  in  the  form 


=‘/ 


where  now  L  is  the  3  x  3  motion  matrix  defined  bv 


K  ’ 

A{n) 


(26) 


T  = 


0  sin  0  0 

^^Aee 

A 


(27) 


We  then  accumulate  observation  samples  in  a  manner  analogous  to  liq.  (19). 


‘7’ 

=‘ii 

L, 

K 


(28) 


Solutions  to  Eq.  (28)  in  the  last-squares  sense  then  provide  a  more  robust  estimate  of  each  scattering 
location  (x„.y„.z„}. 
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5.  SIMULATION 


We  now  develop  an  example  showing  the  utility  and  methodology  used  to  implement  the  enhaneed 
snapshot  imaging  equations.  Consider  the  eone-like  target  geometry  illustrated  in  Figure  8. 


Figure  8.  RV  s im ulat ion  geometry. 


X-band  wideband  seattering  data  having  1  GHz  bandwidth  were  simulated  using  an  1 1 -point,  point- 
seatter  model  (without  shadowing),  with  eaeh  seatterer  having  unity  amplitude  and  (fixed)  random  phase. 
For  illustrative  purposes,  the  seattering  phenomena  ineludes  both  slipping  seatterers  on  the  base  and  sides 
of  the  cone,  and  point  seatterers  on  the  fins,  as  well  as  a  point  seatterer  located  just  below  the  base  and  on 
the  "'nose/'  Slipping  seatterers  arc  particularly  interesting  because  they  exhibit  no  Doppler  with  changes  in 
roll  angle  (e.g.,  azimuthally  symmetric  slots,  ridges,  and  edges).  The  assumed  motion  parameters  and 
associated  6(t),  (t)(t)  motion  used  to  generate  the  data  sequence  E*{f^,t^}-.  where  denotes  the  scattered 
field,  arc  tabulated  in  Figure  9  along  with  plots  of  the  associated  motion  0(t),  (t)(t).  Note  the  rapid  change  in 
(t)(t)  (shown  in  '"dotted  red”  and  clipped  over  the  plot  region)  vs.  the  slow  change  in  0(t)  characterized  by 
f^»fp.  Of  particular  interest  for  this  target  is  the  Nyquist  sampling  limit  defined  by  Eq.  (11).  For 
target  length  --1.5  m  at  near  "‘nose-on”  aspect,  Ae^  ~  0.75  deg. 
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A  typical  range-Doppler  image  used  to  extract  the  sequence  {^n’^nlq  is  also  shown  in  Figure  9. 


Parameter 

Simulation 

K  (deg) 

14.4 

Bp  (deg) 

2.71 

T,  (sec) 

0.63 

Tp  (sec) 

4.0 

Wp  (deg) 

-62 

ae.  (deg) 

158 

T  (sec) 

0.544 

Phi  (0,  2n)  Theta 
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00  . 
oD  0  ^ 

8 


-3-2-10  1  2  3 


Range  Rate  (m/s) 


Figure  9.  Simulation  motion  parameters  6  (t),  (j>(t)  and  typical  range-Doppler  image. 


Figure  10  illustrates  a  difTerent  characterization  of  the  variation  of  [0(t),  (|)(t)]  for  the  assumed  motion 
over  several  roll  cycles,  where  we  plot  0(t)  vs.  (t)(t)  parametrically  vs.  time.  The  “arrows”  on  the  blue  trace 
illustrate  the  direction  change  vs.  time.  The  “phi-euts”  illustrated  represent  typical  theta  values  at  (j)  =  (j)  + 
2*n*7i  jumps.  We  clearly  see  that  for  the  value  of  0p  considered,  for  some  eases,  the  change  A0  between 
A(j)  =  271  Jumps  is  greater  than  Nyquist.  However,  in  the  region  0  =  0,  i.e.,  along  the  upper  portion  of  the 
trajectory,  the  change  delta  theta  is  below  Nyquist  for  the  assumed  0p  =  2.71°.  Thus,  for  example,  the 
processing  region  indicated  on  the  figure  is  well  suited  to  apply  phase-enhanced  3D  imaging.  For  other 
regions,  conventional  snapshot  imaging  techniques  would  be  applicable. 

By  extending  the  above  discussion,  and  using  the  small  precession  expressions  for  O(t),0(t)  defined 
by  Eq.  (8),  one  ean  readily  show  that  when  operating  about  the  region  0  =  0,  the  precession  angle  0p  is 
restricted  by  the  more  general  relationship. 
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Figure  10.  Oft)  vw.  (jXt)  illustrated parauietrlcally  iu  thue. 

I  he  above  diseussion  leads  to  a  methodology  one  might  apply  for  eases  when  the  precession  angle 
0p  is  eonsiderably  larger  than  Nyquist: 

1 .  Isolate  regions  about  6  ~  0  . 

2.  Choose  a  value  (])  =  Seareh  on  t^^  for  (])  =  (j)  +  2jr  according  to  l^  -  Icf  )  -  F  27r\^ 

3.  IVocess  the  associated  sequential  pairs  of  images  for  the  times  C,  ^  according  to  Fq.  (28). 

4.  Sequentially  accumulate  estimates,  by  sliding  t^^  in  steps  At,  keeping  within  the  bounds  of  the 
Nyquist  sampling  criteria. 

I  hc  above  methodology  applied  to  the  simulation  example  eonsidered  leads  to  a  processing  interval 
indicated  as  the  red  portion  of  the  [0,  ([)]  grid  illustrated  in  Figure  10. 
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In  Figure  1 1,  we  illustrate  a  magnifieation  of  the  [9,  0]  sampling  region  eorresponding  to  the  red 
portion  of  Figure  1 0,  indieating  the  0,  0  values  assoeiated  w ith  t^  ((j)  =  0)  and  tq  ((])  =  0  -h  2k).  Clearly,  the 
ehange  A0  over  this  region  is  less  than  Nyquist,  whieh  for  the  aspeet  angle  eonsidered  eorresponds  to 
approximately  0.75®.  Also  shown  in  the  Figure  are  the  (R.  R)  values  extraeted  from  the  range-Doppler 
images  at  time  t  =  t^  as  t^  progresses  over  the  indieated  region. 
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Figure  IL  Image  grid  and  extracted  poles. 

In  Figure  12,  we  illustrate  the  eomposite  of  these  (R,  R)  pairs  for  each  of  the  sliding  intervals  [t^] 
and  [tq].  Note  the  elose  eorrespondenee  of  the  range-Doppler  images  as  eaeh  seatterer  appears  only 
slightly  displaced  due  to  the  change  in  aspect  A0.  Note  also  from  Figure  1 1  that,  in  faet,  A0  changes  with 
time  increment  as  t^  is  incremented  over  the  processing  interval. 

Figure  13  illustrates  the  3D  image  obtained  when  applying  the  phase-enhaneed  snapshot  equations 
to  the  (l^n’  ^  (phase)  J  sequence.  Note  that  because  the  technique  provides  an  independent  estimate  of 
the  actual  3D  location  of  eaeh  seatterer  as  contrasted  to  a  conventional  range-Doppler  image  where  the 
location  of  an  image  peak  must  be  estimated,  the  3D  "‘images”  presented  here  are  actually  estimates  of  the 
3D  locations  of  the  individual  scattering  centers.  The  amplitudes  and  frequency  eharaeteristies  of  eaeh 
seatterer  are  extracted  separately,  as  discussed  in  [1]. 
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Figure  12.  Raiige-Doppler  pairs  (I  Iq  vs.  /  i^). 
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Figure  13.  Phase-enhanced  3D  snapshot  image. 
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Two  facts  arc  immediately  obvious  upon  examination  of  the  results  of  Figure  13. 

The  phase-augmented  snapshot  equations  correctly  image  the  slipping  scatterers. 

The  slipping  scatterers  arc  imaged  on  the  cone  only  over  the  observation  interval — in  this  case  the 
roll  angle  (J)  given  by  28"^  <  <J)  <  92°. 

By  comparison.  Figure  14  illustrates  a  3D  target  image  using  conventional  snapshot  imaging  (i.e.,  no 
phase)  over  a  region  where  A0  >  Nyquist  when  ([)  =  (})  +  27:  as  indicated  on  the  insert.  Over  this  observation 
region,  the  precession  causes  too  large  a  change  in  0  over  one  roll  cycle  so  that  the  extracted  phase  is 
ambiguous,  and  conventional  snapshot  image  techniques  must  be  used.  Note,  however,  that  in  this  case  the 
slipping  scatterers  are  not  properly  located  and  project  to  the  z-axis. 
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Figure  14.  Conventional  snapshot  imaging  (aspect  sample  >  Nyquist). 
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6.  APPLICATION  TO  BISTATICS 


Upon  close  examination  of  the  development  in  Seetion  3.0,  it  beeomes  evident  that  the  phase 
information  contained  in  interrogating  the  target  from  direetions  k,  and  k^  need  not  be  obtained 
sequentially.  The  general  development,  with  limited  modifieation,  holds  if  the  baekseatter  is  obtained 
from  two  sensors  simultaneously,  using  either  monostatie,  bistatie,  or  multistatie  data.  However,  for  the 
ease  of  more  than  a  single  sensor,  the  issue  of  eoherenee,  or  phase  adjustments  to  the  data,  must  be 
addressed.  We  now  develop  the  modifieations  required  for  the  bistatie  ease  for  two  sensors  labeled  "M"  and 
For  eonvenienee,  we  assume  sensor  1  operates  monostatieally,  and  sensor  2  reeeives  the  bistatie  data 
seattered  from  the  incident  field  of  sensor  1. 

The  backscattcred  field  in  the  direction  of  sensor  1  (direction  vector  k ,  )  is  characterized  by 
F^q.  (12),  except  now  we  must  add  an  arbitrary  one-way  propagation  delay,  denoted  as  r,  .  For  this  case 


B; 


-](o{2  r, ) 


I  „  D  ,,  e 


(30) 


In  a  similar  manner,  the  bistatie  Held  in  the  direction  k  ,  is  given  by 
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k,)*  r  , 


(31) 


where  again,  we  assume  for  k,  close  to  k  ^  '  U  =  D  ^  .  Fxpressing  kj  =  ki  ^  A  k  ,  I^^q.  (27)  can  be 
written  in  the  form 
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(32) 


Comparing  (30)  and  (31),  after  extracting  the  differences  in  image  phase  for  the  nth  scatterer,  we 
obtain  the  modified  version  of  Eq.  (20): 


A„  (phase)=|^(AR)+  Ak-r 
where  we  have  defined  the  coherence  parameter  AR  ■ 


(33) 
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AR  -  c  (fj  -  T|)  ’ 


(34) 


and  c  =  3  X  10*^  m/s.  Note,  for  the  bistatic  case,  the  occurrence  of  the  phase  factor  (27r/^)  vs.  (47r/>.)  in 
Bq.(33). 

Clearly,  when  using  Eq.  (32)  to  augment  the  snapshot  equations,  AR  must  be  estimated.  One  ean 
then  either  eorreet  the  data  aeeordingly,  or  use  AR  in  Eq.  (33)  direetly.  Teehniques  for  estimating  AR  are 
deseribed  elsewhere  [13]. 

The  separation  required  between  k^and  k2  to  afFeet  less  than  Nyquist  sampling  is  direetly 
dependent  on  the  range  and  projeeted  length  of  the  target.  Some  typieal  baseline  values,  normalized  to  the 
Nyquist  sampling  interval  A0^  ,  are  illustrated  in  Figure  15. 


Figure  15.  Mciximim  baseline  for  bistatic  collections. 
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7.  SUMMARY 


Three-dimensional  inlerferomelrie  SAR  imaging  of  ground  reflecting  surfaees  has  shown  mueh 
promise  for  synthetie  aperture  radar  proeessing.  A  key  factor  in  achieving  such  3D  images  is  precise 
eontrol  over  i:hc  sampling  grid  (motion)  used  by  the  radar  sensor.  This  report  develops  a  framework  based 
on  3D  snapshot  imaging  of  sparse  angle  sectors  of  ISAR  data,  extending  these  results  to  3D  interferometric 
ISAR  imagirg  of  targets  having  dctciministic,  torque-free  exoatmosphcric  motion  characterized  by  spin 
and  precession. 

Essential  to  the  application  of  these  techniques  to  Eulcr-type  dynamical  motion  is  the  existence  of  a 
sparse-angle  sampling  grid  over  which  the  phase  information  extracted  from  range-Doppler  images  can  be 
properly  utilized.  We  compare  several  types  of  sparse-angle  sampling  grids  to  that  generated  from  Eulcr- 
type  motion  and  develop  an  ISAR-based  approach  to  extracting  the  required  phase  infonnation. 

The  framework  presented  here  is  an  extension  of  the  '‘3D  Snapshot  Imaging'"  technique  developed  in 
1 1 1.  A  new  set  of  independent  snapshot  equations  arc  developed,  from  which  image  points  on  the  target 
can  be  extracted.  Because  the  difTerential  phase  extracted  between  suitable  range-Doppler  images  is 
required  to  be  nonambigiioiis,  the  applicability  of  the  generalized  technique  requires  the  precession  angle 
Op  be  relatively  small  [see  Eq.  (28)|.  However,  the  noncoherent  snapshot  image  equations  remain  valid 
beyond  this  precession  limit,  and  the  unique  sampling  grid  developed  is  still  very  useful  for  3D  imaging  of 
the  noncoherent  snapshot  equation  set.  Techniques  for  removing  this  constraint  by  resolving  the  phase 
ambiguities  for  A0>AO,^  are  under  investigation. 

A  simple  simulation  is  developed  illustrating  the  application  of  the  technique,  and  extensions  to 
bistaties  are  developed.  Applications  of  this  technique  to  field  data,  not  presented  in  this  report,  have 
demonstrated  good  results. 
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